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ALMOST PERIODICITY IN TIME OF SOLUTIONS OF THE KDV EQUATION 


ILIA BINDER, DAVID DAMANIK, MICHAEL GOLDSTEIN, AND MILIVOJE LUKIC 


Abstract. We study the Cauchy problem for the KdV equation dtu — ()udxU+d^u = 0 with almost periodic 
initial data 11(3;, 0 ) = V{x). We consider initial data V, for which the associated Schrodinger operator is 
absolutely continuous and has a spectrum that is not too thin in a sense we specify, and show the existence, 
uniqueness, and almost periodicity in time of solutions. This establishes a conjecture of Percy Deift for 
this class of initial data. The result is shown to apply to all small analytic quasiperiodic initial data with 
Diophantine frequency vector. 
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1. Introduction 

This paper is devoted to the Cauchy problem for the KdV equation, 

( 1 . 1 ) dtu — 6udxU + d^u = 0, 

(1.2) u{x,Q) = V{x) 

with initial data V : R — R. It is assumed that u(x, t) : R^ —R and solutions are considered in the classical 
sense: u is three times differentiable in x, once differentiable in t, and obeys (ED for each (x, t) G R^. 

The KdV equation is named after Korteweg-de Vries |KdV] . who proposed it as a model for the propa¬ 
gation of shallow water waves in a canal. It has received a lot of attention after the discovery in 1960s of 
infinitely many conserved quantities and the inverse scattering transform by Gardner-Greene-Kruskal-Miura 
[GGKMl IGKM] and the Lax pair formalism |Lax] which showed isospectrality of the associated Schrodinger 
operators —d^ + Due to these properties, the KdV equation is often described as completely inte- 

grable. This formalism was first implemented for sufficiently fast decaying initial data U, where the conserved 
quantities are integrals of polynomial expressions in the solution and its derivatives, the lowest one being 
the (square of the) norm, Jg\u{x,t)\'^dx. 
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In the 1970s, this framework was successfully extended to solutions of the KdV equation on a torus, 
It : T X R — 7> R, where the conserved quantities, of course, take the form of integrals on the torus, e.g. 
Jj\u{x,t)\'^dx. It was shown, through the contributions of many authors [Gal IFZl iDull INo[ IMcKvMl IDMNl 
IFMcLl IMcKTI IDu2) , that in this setting, the KdV equation is a completely integrable Hamiltonian system 
with action-angle variables; for a book treatment, see [KPj . In particular, the existence of action-angle 
variables implies that solutions are almost periodic in t. Of course, solutions on T x R can be viewed as 
spatially periodic solutions on R x R, which motivates the following question. 

Deift [De] posed an open problem whether, for almost periodic initial data V, the solution of dm, doj is 
almost periodic as a function of t. The analysis of the KdV Cauchy problem with almost periodic initial data 
presents significant new obstacles. The conserved quantities in this setting would take the form of spatial 
averages, e.g. 


lim —— 
L—>oo 2T 


\u{x^ t)\'^dx, 


and it is not known how to use such conserved quantities to prove even short time existence of a solution. 
Moreover, isospectrality of Schrbdinger operators is too weak a condition to restrict the phase space of our 
solution to a torus, a first prerequisite for true integrability. In addition, unlike Schrbdinger operators with 
fast-decaying or periodic potentials, almost periodic Schrbdinger operators have diverse spectral properties. 

In this paper, we answer Deift’s question for the case of small quasiperiodic analytic initial data with 
Diophantine frequency. Let us first define the appropriate class of quasiperiodic functions and then state the 
theorem. 


Definition. Let e > 0, 0 < Kq < 1, and w G for some v G N. We say that V G J’(a;, e, Kq) ifV : R —>■ R 
is of the form 

(1.3) V{x) = , 


such that 

(1.4) |c(n)| < eexp(—Ko|n|), 

We make J’(w,e, kq) a metric space with the metric inherited from L°°(M). 

Theorem 1. Let uj G obey the Diophantine condition 

(1.5) |nw| > ao|u|“^“, u G \ {0} 
for some 

(1.6) 0 < oo < 1, V < 6o < oo. 

There exists £o{ao,bo, kq) > 0 such that, if e < Eq and V G CP(w,e, kq), then: 

(i) there is a unique solution u of (II.II) . (11.21) such that 

(1.7) u,a,,,uG L“(Mx [-T,r]), VT<oo. 

( a ) for every t G R, u{-,t) is quasiperiodic; in fact, u{-,t) G CP(a;, \/4e, ko/4) 

(Hi) u is almost periodic in t, in the following sense: there is a compact (finite or infinite dimensional) 
torus a continuous map 

M : ^ T(w,\/4i, ko/4), 

a base point a G T'^, and a direction vector ( G R^^ such that u{-,t) = M(q; -|- (t). 

Throughout this text, T = R/27rZ. 

In (11.41) . (II. Sp . and elsewhere in this work, |n| := ^ ^ ■ 
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Remark. On Kq), the L°°-norm is equivalent with the norm 


11 ^-^ 11 . 



|2g2|n|r 


1/2 


for any r < kq, and with the Sobolev norm inherited from Vk^’°°(]R) for any fc G N. Part (Hi) of the previous 
theorem therefore implies that besides u, derivatives ofu are also almost periodic in t, and so is each Fourier 
coefficient c(n,f) ofu{x,t). 


The initial value problem (HID, (HID for quasiperiodic initial data was previously approached by Tsug- 
awa [Ts| and Damanik-Goldstein [DG2) , building upon work of Kenig-Ponce-Vega [KPV] and Bourgain [Bo] . 
[Ts] proved local existence and uniqueness for initial data of the form dESD with Diophantine u) and 
c(n) decaying at a sufficiently fast polynomial rate, and |DG2) proved local existence and uniqueness for 
V G J’(w, e, Kq) for all uj and global existence and uniqueness for Diophantine uj. However, in all those 
results, uniqueness is proved within a class of solutions such that u(-, t) is quasiperiodic for all t, with similar 
decay conditions on Fourier coefficients as those assumed for the initial data (polynomial in [Ts] and expo¬ 
nential in [DG2) 1. Our Theorem [T] improves on this by proving uniqueness within a larger class of solutions 
assuming a priori only a boundedness condition (ED and not quasiperiodicity. Moreover, while the class of 
initial data in Theorem [T] is the same as that considered in |DG2) . this paper provides an alternative proof of 
the results from |DG2) for Diophantine uj. The alternative proof avoids the detailed combinatorial analysis 
in [DG2) and relies instead on the inverse spectral theory of Schrodinger operators and fundamental work of 
Rybkin |Ry| , coupled with [DGL2] . Finally, Theorem [T] proves almost periodicity of the solution in t, which 
was beyond the reach of the approach in [Tsi IDG2] . 

In fact, we are able to prove a more general result, and prove existence, uniqueness, and almost periodicity 
in t whenever V is almost periodic and the spectrum of the associated Schrodinger operator —df -|- H is 
absolutely continuous and not too thin, in a sense to be made rigorous below. Stating that result requires 
introducing additional prerequisites. 

For W : R —^ R, consider the Schrodinger operator 


(1.8) [Hwy]ix) = -y"{x) FW{x)y{x), a; G R. 

In the cases of interest in this paper, W will be a bounded function, so (11.8(1 defines an unbounded self-adjoint 
operator on L^(R) with the domain 

D{Hw) = D{Ho) = {y : R^ R I y G AGioc(R),y' G AGioc(R),y" G L‘^{R)}. 


The spectrum a{Hw) is the set of 2 : G C for which Hw — z does not have a bounded inverse map {Hw — z)~^ 
L^(R) —>■ D{F[w). Since FIw is self-adjoint, a{Hw) C R. 

For if G T^(R), the spectral measure dp.^ is the unique measure on R with the property that 


(V', {Hw 



- dpHx), 

X — z 


Vz G C \ R. 


Using the Lebesgue decomposition of = dy^,ac + d/x^_sing with respect to Lebesgue measure, the abso¬ 
lutely continuous spectrum of Hw can be defined as the smallest common topological support of absolutely 
continuous parts of all spectral measures. 


<Xb.c{Hw)= U suppdy^^ac- 

V>GL^(R) 

Glearly, suppd/x^^ac C suppd/i^ C a{Hw), so o'ac(Rvu) C a{Hw). We will later be interested in cases in 
which as,c{Hw) = cr{Hw). 

The spectrum S = a{Hw) is closed and bounded from below but not from above, so it can be written in 
the form 

^=[A,oo)\ 

jeJ 
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where S = inf S and {Ej , ) are the maximal open intervals in R \ called gaps. We will sometimes also 

work with the closure of S in the Riemann sphere, denoted S = S U {oo}. 

Denote by gc\si^) the potential theoretic Green’s function for the domain C \ S with a logarithmic pole 
a,t E— 1. The set S is regular if 


lim oc\s( 2 ) = 0 , \/xGS. 
zec\s 


(1.9) 

Denote by Cj S {E~,E^) the critical points of gc\s in the gaps. The Parreau-Widom condition states that 


( 1 . 10 ) 


J^ScXsicj) < oo. 
jeJ 


Regularity and the Parreau-Widom condition were first used in spectral theory in groundbreaking work 
of Sodin-Yuditskii [SYll ISY2) . The results in those papers are stated in terms of homogeneity in the sense 
of Carleson, that is, the existence of r > 0 such that 


( 1 . 11 ) 


ISTl [a;o - e,a;o + e]| > re, Vxq G S', VeG(0,1], 


and finite gap length. However, it is then remarked that those conditions imply regularity and the Parreau- 
Widom condition, and the remainders of their proofs rely only on that. 

We will need additional conditions on the spectrum. We denote 


for j € J and 




r?,,, = dist((i?-, £;+),(£;,-, £;,+)) 

gjfi = d\si{{E~ ,E^),E) 


3 ’ 3 

for j,l G J (notationally, we assume here that our abstract index set J does not contain 0 as an element). 
We also denote 

(1.12) Q = (r?,-o + 7,r/"n(i + ;r'^ 


leJ 

1¥^3 


We assume that S satisfies a set of conditions of Craig |Cr) : 

7 I /2 

(1.13) ^ 7 j < 00 , sup 7 j( 7 j < 00 , sup ——Cj < 00 , 


3^J 


3&J 


3^J Vjfi 


1/2 1/2 

supY^ ————Cj < 00 . 
3^jf^j V3,l 
1¥^3 


Craig introduced these conditions to ensure that a vector field, associated with the translation flow, is 
Lipshitz on the isospectral torus T>{S) = with the metric given by 


(1.14) 


11/5 - (^IId(s) = sup7j/^||(/)j - (fjWj. 


3&J 


In order to ensure that a second vector field, associated with the KdV flow, is Lipshitz, we need to 
strengthen these conditions as follows: 

(1.15) 




7/ <00, 


sup 7 


i - —Cj < 00 , sup y — - (1 + r]j,o)Cj < 00 for a G 1 > . 

V 3,0 \ ^1’' / J 

1^3 

Note that while these Craig-type conditions look complicated, they are quite weak. For instance, it is 
straightforward to see that they are satisfied for the small quasiperiodic initial data of Theorem [T] Also note 
that (11.151) implies (11.131) . 
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We also wish to emphasize that this form of Craig-type conditions is sufficient, but not necessary. In 
particular, if one changed the metric (I1.14p . the form of the conditions (11.131) . (jl.l5l) would change. At 
this stage in the theory, the fine structure of gap sizes and distances is not well understood for general 
almost periodic operators, so it is not clear that varying the metric would broaden the scope of the following 
theorem. We therefore choose, for clarity of exposition, to keep the metric (11.1411 used by Craig. 

To ensure continuity of scalar fields associated with certain trace formulas, we will also assume that 

(1.16) ^(1+ 7?|o)7j < oo. 

Theorem 2. Let the initial data 1^ : K —>■ R 6e uniformly almost periodic. Denote S = a{Hy) and assume 
that 

(i) S = Crac(i?V'); 

(a) S is a regular Parreau-Widom set; 

(Hi) S obeys the Craig-type conditions (11.1511 . (ll.lfiH .- 
Then: 

(i) there is a unigue solution u of (ll.l|l . (11.211 which obeys the boundedness condition (IE3; 

(a) the solution u is almost periodic in t, in the following sense: there is a continuous map 

M : ^ W‘‘’“(R), 

a base point a € and a direction vector ( G R'^ such that u{-,t) = At(a -I- ft); 

(Hi) for each t S R, the function u{-,t) is uniformly almost periodic with freguency module egual to the 
frequency module ofV. 

We recall that the frequency module of an almost periodic function V is the Z-module generated by the 
set of all A G R such that 

lim ^ f V{x)e-^^'-^dx^0. 

L —^oo ZiL J _ 

The condition (EH in Theorems [T] and [5] is assumed only in order to use the results of Rybkin |Ry| ; as 
discussed in Remark 2 on page 9 of |Ry| , this condition can be relaxed. 

A construction of Cohen-Kappeler |CK] shows nonuniqueness of the Cauchy problem (11.111 . (11.21) even for 
y = 0 on some domains of the form {{x, t) G R^|0 < t < h{x)} with h a strictly positive, increasing function 
of X. 

An affirmative answer to Deift’s problem was already known in some cases. A class of “algebro-geometric” 
solutions of the KdV equation was constructed in the 1970s, resulting from the same theory which solved 
the periodic case. These solutions correspond to almost periodic Schrodinger operators with finitely many 
gaps in the spectrum, so they are often called finite zone or finite gap solutions. Due to having finitely many 
gaps, their frequency module under any linear flow is finitely generated, so they are quasiperiodic in both x 
and t. We emphasize, however, the difference between these finite gap solutions and our quasiperiodic initial 
data of Theorem (TJ which generically have infinitely many gaps in the spectrum. 

A result of Egorova |Eg| goes beyond the finite gap case to construct almost periodic solutions with 
gap sizes obeying a superexponential decay condition in a parameter resembling the period of periodic 
approximants. This result applies to a class of limit-periodic potentials previously studied by Chulaevsky [Ch] 
and Pastur-Tkachenko [PTll IPT2) , which can be approximated superexponentially in the period by periodic 
potentials. Namely, V is assumed to have a„-periodic potentials 14, such that for any 5 > 0, 

lim sup||y - Vn\\w^,mx,x-\-i)) = 0. 

xeK 

Due to the superexponential rate of approximation, this analysis can obtain convergence from various esti¬ 
mates exponential in the period. This approach would not apply in our setting. 

While Theorem [1] is direct in the sense that it derives the desired conclusion from explicit assumptions 
about the initial datum, there is a sophisticated inverse spectral theory working in the background which 
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makes this result possible. The inverse spectral theory depends on fundamental developments by Craig [Crj . 
Sodin-Yuditskii [SYlj . Remling [Relj . centered on the notion of reflectionless Schrodinger operators. Another 
fundamental development is due to Rybkin |Ry| , who studied the time evolution of solutions of dm in terms 
of the Weyl M-matrices of the associated family of Schrodinger operators; since the M-matrix determines a 
Schrodinger operator uniquely, Rybkin’s approach provides much more information about the solution than 
what can be concluded merely from unitary equivalence of operators. 

We wish to emphasize that Theorem [5] is not itself a small coupling result, even though its corollary. 
Theorem [3 is. It is natural to ask to what extent these properties persist as the coupling constant e is 
increased. In this regard, a particularly well understood example is the almost Mathieu operator, a discrete 
Schrodinger operator on given by {H\u)n = + m„+i + 2A cos(27r(na; + 9))um where A > 0, 

a; G M \ Q, 0 G K. The almost Mathieu operator goes through a transition from absolutely continuous 
to singular spectrum as the coupling constant A goes through 1. Moreover, the Lebesgue measure of the 
spectrum is 4 |A — 1|, so A = 1 is also a critical point in that sense. Many authors have contributed to these 
results; see the recent survey papers |Da[ IJiMa) and references therein. Recently, it was shown in |DGSV] 
that the spectrum of the almost Mathieu operator is homogeneous whenever A 1 (under a Diophantine 
condition for w). This is encouraging evidence that for analytic quasiperiodic operators, the assumptions of 
Theorem [2] may hold as we increase the coupling constant, up to a critical value where there is a transition 
in the spectral type. However, this is currently in the domain of speculation. 

More specifically, our work serves as motivation for one to work out a continuum analogue of Avila’s global 
theory for discrete one-frequency analytic quasi-periodic Schrodinger operators |Av IE V IE especially 
of his results in the subcritical regime. That is, one should classify cocycle behavior in the continuum case in 
the same way as Avila does in the discrete case, namely as subcritical/critical/supercritical, and prove that 
subcritical cocycle behavior leads to almost reducibility and absolutely continuous spectrum. The ultimate 
goal would be to show that our results apply throughout the subcritical regime. Of course, one would also 
have to establish the necessary structural properties for the spectrum as a set, and here the result from 
[DGSV] suggests that these properties may indeed be typical in the subcritical regime. 

While the previous paragraph describes ways in which the known results in the discrete case may inform 
a study of the continuum case, one can also turn this around and attempt to advance the theory in the 
discrete case by implementing ideas that have proved to be useful in the continuum case. For example, given 
the results obtained in the present paper, the following natural question arises: considering the Toda lattice, 

d(n, t) = a{n, t){b(n -|- 1, t) — b(n, t)), 

(1.17) 

b{n, t) = 2(a^(n, t) — a^{n — 1, t)), 

so that the initial data (a(-, 0), &(•, 0)) are almost periodic and the associated Jacobi matrix is reflectionless 
and its spectrum is a regular Parreau-Widom set and satisfies suitable Graig-type conditions, is it then true 
that there is a global solution that is almost periodic in time? Based on the close analogy between the KdV 
equation and the Toda lattice, as well as our Theorem [2l we expect an affirmative answer to this question. 
Many of the key ingredients are already in place; for example, Remling has shown the reflectionlessness of 
almost periodic Jacobi matrices with fully supported absolutely continuous spectrum [Ri2] and the invariance 
of reflectionlessness under the Toda flow [Re3] , and Sodin and Yuditskii’s work exists in the discrete setting 
as well |SY2] . 

An attractive special case is of course given by almost Mathieu initial data, 

a(n, 0) = 1, b(n, 0) = 2Acos(27r(nw -|- 0)). 

This special case is particularly interesting because, as indicated above, the spectral analysis of the initial data 
is almost complete. In particular, the phase transition at A = 1 is understood on a variety of levels: Lebesgue 
measure of the spectrum, quantum transport, spectral type, behavior of (generalized) eigenfunctions, and 
more concretely the behavior of the underlying cocycles. On the one hand, this suggests that one may 
hope that by working out the discrete analogue of our work, one can establish the desired result on almost 
periodicity in time throughout the subcritical region (i.e., for A < 1). On the other hand, this case exemplifies 
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the need for significant additional ideas to penetrate the supercritical region (A > 1). Of course, it is known 
that in the supercritical region, the initial data will not be reflectionless. However, the structure of the 
spectrum as a set is the same as in the subcritical region; indeed, the spectra corresponding to coupling 
constants A and A“^ are related via linear scaling due to Aubry duality. This shows clearly that isospectrality 
alone is insufficient to get a good understanding of the evolution (even with additionally imposed Craig-type 
estimates). At the present time, it is unclear how to define the appropriate “isospectral torus” in the 
supercritical regime, or more generally when reflectionlessness fails. The critical regime (in the sense of 
Avila, corresponding to the case of critical coupling, A = 1, in the case of the almost Mathieu operator) 
poses additional challenges due to the critical nature of the associated cocycles and the degenerate structure 
of the spectrum. 

While the critical regime is atypical in the class of discrete Schrodinger operators with analytic one- 
frequency quasi-periodic potentials, the features of the critical regime, such as spectra of zero Lebesgue 
measure and purely singular continuous spectral measures, are actually generic in the continuous category 
(i.e., varying the sampling function in (7° while keeping the base transformation fixed). For example, it was 
shown by Avila [Avl] in the discrete case and then, based on his ideas, by Damanik-Fillman-Lukic [DFL] in 
the continuum case that the spectrum of a generic bounded continuous limit-periodic operator has Lebesgue 
measure 0. A result of Avila-Damanik [AD] shows that for discrete quasiperiodic Schrodinger operators, 
singular spectrum (and hence the absence of reflectionlessness!) is generic, and this suggests that the same 
is true for continuum quasiperiodic Schrodinger operatorsQ (It is conjectured that zero-measure spectrum 
is generic in this setting as well; however, this is not yet known, see [ADZ] for a partial result.) Finally, for 
discrete Schrodinger operators, the absence of point spectrum was also shown to be generic in the continuous 
category by Boshernitzan-Damanik m- These results suggest that a complete answer to Deift’s problem 
will require the development of robust new tools in inverse spectral theory. 

In any event, the spectral analysis of almost periodic Schrodinger operators and Jacobi matrices is in 
general very complicated and just having a Lax pair is not enough to integrate the evolution equation (i.e., 
the KdV equation and the Toda lattice, respectively) 

As a final remark in this introduction, let us point out that we obtain the almost periodicity of the solution 
u as a function on Namely, instead of writing = M(a -I- (t) as in Theorem [21 we can also write 

u{x, t) = F{a + -I- Qt) with a continuous function F : ^ R. This follows from our proof of the identity 

u{-,t) = M(a -I- Qt) (along with a suitable composition with the Dubrovin flow); compare Section 0 

Sections|2]and|3]prepare the stage for the uniqueness part of Theorem|2] Section|2]starts with the necessary 
definitions and reviews and expands previous results on the evolution of Dirichlet data of a reflectionless 
operator under translation, and Section [3] develops a similar description for the evolution of Dirichlet data 
under the KdV flow. Section 2] is devoted to the existence part of Theorem |21 and Section |S] to the almost 
periodicity of the KdV flow. Section [5] also contains the proof of Theorem |21 based on all the previous 
material. Finally, Section [6] uses Theorem [2] to prove Theorem [TJ 

2. The Translation Flow and Non-Pausing of Dirichlet Eigenvalues 

In this section, we discuss the translation flow applied to a reflectionless Schrodinger operator with spec¬ 
trum S obeying (jl.ljp . In particular, this flow can be expressed in terms of the Dirichlet data of the operator. 
We begin by introducing the notion of reflectionlessness and reviewing results of Craig [Crj . after which we 
investigate the behavior of Dirichlet data at gap edges to give a new characterization of points where a 
Dirichlet eigenvalue goes through a gap edge and, in particular, to explain why Dirichlet eigenvalues don’t 
pause at gap edges. 

For z G C \ a{Hw), the second order differential equation 

-y" + Wy = zy 


^We expect that it is not too hard to work out a continuum analog of this particular result from [ADj. Given the obvious 
relevance to the scope of our work in the context of the KdV equation, it would be nice to see this extension be worked out 
explicitly. 
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has nontrivial solutions tj}±{x;z), called Weyl solutions, such that ±oo), dx). In terms of 

the Weyl solutions, one can express the half-line m-functions associated with the half-line restrictions of Hw 
to [x, ±oo) with a Dirichlet boundary condition at x by 


m±{x; z) 


'4’±{x;z) 

'4’±{x;z)' 


For each x, these are meromorphic functions of z € C \ a{Hw)- 
The diagonal Green’s function associated with Hw is given by 


G{x, x; z) 


_ i’-jx; z)il;+{x; z) _ 

tp'_ix; z)ip+{x; z) - ijj'^{x; z)ijj-{x; z) 


1 

TO_(a;; z) — m+{x; z) ’ 


and it is an analytic function of z G C \ a{Hw) for each x. 

The function G{x, x; z) has nontangential limits G(x, x]y + iO) for Lebesgue-a.e. y G M, and Hw is called 
reflectionless if 

3?G(x, x-,y + iO) = 0 for Lebesgue-a.e. y G a{Hw)- 

It is known [SYI] that this definition is independent of the choice of x. We denote the set of reflectionless 
operators with spectrum S by 3?(5'). 

Recall the set T){S) = with the metric (11.141) . We introduce variables on T){S) given by 


(2.1) 

H = Ej + {E+ - E. ) cos^{ipj/2) 


f -1- ifj € (0, tt) -I- 27rZ 

(2.2) 

CTj = < — ifj G (—TT, 0) -1- 27rZ 


[o G ttZ 


On T>{S), we introduce a scalar field Qi, 


Qii^)=E + J2iE-+E+-2y,) 

jeJ 

and a vector field iL, with components 

The tangent space of a point on 'D{S) will be equipped with the norm 

(2.3) ||n|| = supyj^Vjl 

j&J 

and vector fields on ©(S') will be equipped with the sup-norm obtained from (12.31) . 

The motivation for this setup comes from the Dirichlet data of a reflectionless operator on S. Let W be 
reflectionless and a{Hw) = S. Following the definition of Craig [Cr], the Dirichlet data Hj{x) are defined as 
zeros of the diagonal Green’s function in the gaps. More precisely, using the observation that G{x, x; z) is 
strictly increasing for z G {E~ ,E^), Craig defines 

(z€(E~,E+) G(x,x;z)=0 

(2.4) y,j[x) = <E^ G{x,x;z)<Oioial\ze{E~,E^) 

[ E~ G{x, x;z)>0 for all z G {E~, E ^) 

If y.j{x) G {E~,E^), the sign aj{x) G {-I-,—} is also defined. It is determined in terms of which of the 
half-line m-functions m+, m_ has a pole at namely, so that nia-j has a pole at y,j. This definition 
immediately implies that for any choice of z G {E~,E^) and a G {-)-, — }, 

(2.5) {a; G M I fJ,j{x) = z,aj(x) = cr} = {x G K | ipaix] z) = 0} 
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where i^±{x; z) are the Weyl solutions at the point z G {E~,E^). 

One then also defines angular coordinates (pj{x) G T implicitly by (12.11) . (12.21) . 

We therefore have, corresponding to the potential W, a trajectory (f{x) on D(5'). In particular, if we fix 
X = 0, we obtain a map B : IR(S') —^ 21(5') given by 


(2.6) B{W) = (/7(0). 

Craig [Cr] showed the following: 

(i) As a function of x G M, ^^(x) is a continuous trajectory on D{S). 

(ii) The potential can be recovered from this trajectory by the trace formula 

W{x) = Qi{ip{x)). 

(ill) The diagonal Green’s function of is given by the product formula 

1 


(2.7) 


G(x, x; z) = 


— /GJ 


n 


(aiz(x) - z )^ 
{E- - z){E+ - z) 


with the branch of square root chosen so that G(x, x; z) is analytic in z G C \ S' and G(x, x; z) >0 for 
z < G. In particular. 


( 2 . 8 ) 


5^G(x,x;z)|^^^.(^) = 


n 


(W - ^^3? 


\{E- - Hj)iE+ - ^j) f-J- {El - fij){E+ - fj.j) 


(iv) dt is a Lipshitz vector field on 2I(S). 

(v) Fix j G J. At any point x G M at which Hj{x) G {E~ ,E^), the function ipj{x) is differentiable and 

(2.9) dx^Pjix) = 'iij{(p{x)). 


(vi) If 31{S) is equipped with the topology of uniform convergence on compacts, the map B is continuous. 

We wish to clarify the situation at the gap edges, with respect to (v) above. We will show that the 
characterization (12.51) has an analog for gap edges, although there are no Weyl solutions at a gap edge. In 
particular, this will show that Hj{x) does not pause at a gap edge and that <Pj{x) is differentiable and (12.91) 
holds at all x. 


Proposition 2.1. Assume that S obeys ()1.13|) and W G 1R(S). Fix a gap {Ej ,E^) and let E G {E^ 

There is a nontrivial solution ip of Ew'ip = Eip which is the limit of Weyl solutions from the gap, in the 
following sense: there are normalizing constants c±{z) G C \ {0} for z G [EJ,E'^) such that 

lim c±{z)ip±{x; z) = ip{x) 

-&(E-,E+) 

z^E 

uniformly on compacts in x. For this solution ip, we have 

(2.10) {x G K I ixj{x) = E} = {x G M I ip{x) = 0}. 

In particular, this set is discrete. 

As pointed out in [Crj . knowing continuity of g:j{x) and knowing that (12.91) holds away from a discrete 
set in X implies that (12.9|) holds for all x. Thus, Dirichlet data of W are trajectories of the flow ik, solving 

( 2 . 11 ) d,,ip{x) = ^{ip{x)). 


Conversely, since is Lipschitz, a solution of (12.111) is uniquely determined by its initial data :/5(0), and then 
so is W. Therefore, the map B given by dm is injective. 
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Proof of Prop. \2.1\ . Separating the term I = j in the product formula (lO) , we write 

- zf 1 


G(x., x\ z) = 


1 


n 


{fll{x) - z)2 


{E+ -z)(E--z) 2^E-Z [Ef - z)iE+ - z) 


The product over I ^ j has a finite limit as z ^ E since the l-th term is bounded by 1 + 'ji/pij. We can 
therefore conclude that 


( 2 . 12 ) 


lim \G{x,x; z)\ 
.6(i5-.£+) 

z^E 


0 = E 

+00 fij (x) ^ E 


We normalize the Weyl solutions for z S (Ej , E '^) by making them real-valued and by requiring 

IV'±(0;2:)|^ + |9x'i/'±(0;2:)P = 1. 

We introduce the Priifer phase 9±{x-, z) G R/ttZ by 

dx'4’± cos0± — ijj± sin0± = 0. 

Using compactness, select a sequence {E~,E^) B Zk ^ E so that the Priifer phases at x = 0 stabilize, i.e., 
the limits 

lim 6'±(0;Zfe) 

k—¥CC) 

exist. By stability of solutions of ODEs, we can then conclude that i/'±(x; Zk) converges in the G^ sense on 
compacts in x to solutions ip± of the equation 

Hw'^± = E'4}±. 


In particular, then, 

(2.13) Yivn 9±{x]Zk) = d±{x) 

k—¥oc 


with 9±{x) G R/ttZ determined by 

dx'4’± cos 9± — ip± sin0± = 0. 
The diagonal Green’s function can be written as 

cos 9-{x;z )cos 9+{x;z) 


G{x, x; z) = 


sin(0_(x; z) — 0+(x; z)) 


so, by (I2.12L we have 


(2.14) 


lim 


.e(E-,E+) 

z^E 


cos 9- (x; z) cos 9+{x]z) 
sin(6*_(x; z) — 9+{x; z)) 


0 = E 

+00 /ij(x) E 


Since w G {E~ , E^ ) is above the bottom of the essential spectrum, tp- (x; w) has infinitely many zeros, so we 
can pick x such that p,j{x) = w G {E~ ,E^). Then, from (j2.14l) and (12.131) we conclude that 0-(x) = 9+{x) 
for that value of x; but then that is true for all x G R, i.e. 1 /)+ = ctp- for some c G {—1,+1}. By a 
renormalization we can assume that V’- = 4’+ ^md denote that solution by ip from now on. 

Now assume that x is such that p,j{x) = E. From (12.141) and (I2.13|) . we conclude that cos0(x) = 0, that 
is, '0(x) = 0. We have therefore shown that 

H := {x G R I /ij(x) = E} C {x G R I V'(x) = 0} =: D. 

In particular, we conclude that D is discrete, since D is. 

Since ipj is continuous, (12.9p holds away from a discrete set, and > 0, we conclude that Lpj is a strictly 
increasing function of x. Therefore, D strictly interlaces the set of zeros of ip-{x',w) for w G {E~,E^). 
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We will now show that D strictly interlaces the same set, using the renormalized oscillation theory of 
Gesztesy-Simon-Teschl |GST] . By Theorem 1.3 of |GST) . if z, w are such that (min(z, w), max(z, w)) (13 = 
0, then the Wronskian of and il)+{x\z) is nonzero for all x. This is equivalent to saying that 

9-{x-,w) — 0+(x;z) ^ ttZ for all x. 

Thus, if we lift 9± to continuous R-valued functions, there is an n £ Z such that for all x, 

mr < 6*_(x; w) — 0+{x] z) < {n + l)7r. 

By continuity in z £ {E~ ,E^), n is independent of z. Taking z = Zk and taking the limit Zk ^ E we obtain 

nn < 0_(x; w) — 9{x) < {n + l)7r 

However, by Prop. 3.4 of |GST] . the function 0_(x; w) — 6{x) cannot have a local extremum where its value 
is in ttZ, so the previous estimates strengthen to 

nn < 0-{x] w) — 6{x) < (n + l)7r 

which means that the zeros of ijj-ix; w) and of strictly interlace. 

Knowing that D G D and that each of those sets strictly interlaces the set of zeros of 'tp-{x;w) proves 
that D = D, since i/’-(x; w) has zeros on any half-line in R. 

Finally, since (12.101) determines ij) uniquely (up to a multiplicative factor), we conclude that the limit '0 
is independent of subsequence Zk ^ E. □ 


3. A Dubrovin-Type Formula for the KdV Flow 


In this section, we derive a formula for the time evolution of Dirichlet data for the KdV flow. This relies 
on ideas and results of Graig |Gr) . discussed in the previous section, and of Rybkin |Ry| , who derived the 
time evolution of the Weyl M-matrix for the KdV flow. 

We will consider a solution u{x,t) : —>■ K of the initial value problem (II.1|) . (11.21) . As it has been 

mentioned above, by the Lax pair formalism, the spectrum S = ()) is independent of t. A much more 

recent result of Rybkin |Ry| shows that the reflectionless property is preserved along the KdV flow. More 
precisely, |Ry| assumes that u{x,t) is a solution of (11.11) such that u, dxU, dtu £ L°°(R x [0,T]) for some 
T > 0, and proves the following: 

(i) The Weyl M-matrix is defined as 


.m-f- 


M = 


m_ —rrijf. 

1 m_ +m+ 
^ 2 m_ —m_|. 


1 \ 

2 m_ —mjf. \ 

] 

m_—7714. / 


where m±(x, t; z) are the Weyl m-functions corresponding to m(-, t) on [x, ± 00 ). The M-matrix evolves 

by 


(3.1) 


where P(x,t; z) is given by 


P = 


dtM = PM -H MP^ 


dxU 2{u — z){u + 2z) — d^i 


2{u + 2z) 


-dxU 


This time evolution is valid for all z £ C\a{Hv)- It is also valid in the sense of the nontangential limit 
M(x, t;z + *0) for those z £ a{Hv) for which that nontangential limit exists for t = 0. 

(ii) For Lebesgue-a.e. z £ <Tac{Hv), the value of the reflection coefficient 


m_|_ — TO_ 


is independent of t £ [0,T]. 


m_|_ — TO_ 
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The condition that the reflection coefficient be equal to 0 for a.e. z € S is equivalent to reflectionlessness 
of the potential; see |SY1) . In particular, Rybkin’s results imply that if Y G 31(5') and 5 = (JadHy), then 
G lk(5) for all t G [0,r]. 

Since is reflectionless for all t, we can introduce the corresponding Dirichlet data ip{x,t) as in the 

previous section. 

Recall from the previous section the torus T>(5) and, on it, the scalar field Qi and vector field if. We 
now introduce the vector field S on T)(S), with components 

— ~2((5i + 

Proposition 3.1. If S obeys (11.151) . then 'E. is a Lipshitz vector field on 1){S). 


Proof. We begin by noting that each component Ej is bounded as 

||s,iu<4(||giiu +21^1+27?,- 0 )^ 

Next, each component Ej has partial derivatives given by 


d^j / Qi + 2fj,j 


difk V Mfc - 


-2 ^'^■7fcSin</5fe, kdj 


1 / 1 


difj -fij ' E+ - Mfc - p.j 


E 




from which it is easy to estimate 


[I] 


d^Pk 

-H 

dE, 

if 

dipj 



Vj.k 


E 


Ik 


kdj 

(llQilloo + 2|S| + 2r]jp) Cj"fj 


and then 




3<SJ 


<sup7j/"E 




kGJ 


\\Pk-‘fk\\ 


difk 


< ||(/j-(^||sup^7j^Sfe^^^ 


dipk 


It remains to show that the sup is finite. Splitting off the k = j term from the sum, the remainder can be 
estimated by 

1/2 1/2 / llQilloo + 2|S| + 2r7j 0 , n\ ri ^ r* i , 1/2^^ 1/2 1 ^ 

sup 27 / 7 ^' ---^ + 2 Q^supCn -(l + ^j-,o)C',+ 7 / < 

fe/j k 'b’fc / \ '/>fc / 


and the k = j term can be estimated by 


k¥=j 


k^j 


sup 

j 


E' 


Ik 


IjAm.k + Ik) 


(llQilloo + 2|^| + 2r]j^Q) Cj'yj < 00 


from the conditions (ll.lbp . 


□ 


00 
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Proposition 3.2. Let u{x,t) be a solution of (11.11) . (11.21) sueh that u, d^u, dtu G x [0,T]) for some 

T > 0. Assume that V G 3^(5') and S = (Tac(^^y)- The function (p{x,t) is jointly continuous in {x,t); it is 
differentiable in t and obeys 

(3.2) dt(pj{x,t) = Ej{ip{x,t)). 

As a first step, we prove a version of this proposition away from gap edges. 

Lemma 3.3. Under the assumptions of Prov. \3.‘A at any {x,t) such that ipj{x,t) ^ ttZ, ipj(x,t) is differen¬ 
tiable in t and obeys dSJ). 


Proof. The M-matrix is analytic in z G C \ 5". It has a removable singularity at a Dirichlet eigenvalue, where 
it is equal to 


(3.3) 


M (x, t; Hj{x, t)) = 


_ f-aj{x,t)m_^.(^a:,t)iUjix,t)) -l<7j{x,t) 


-^aj{x,t) 


0 


Note that the bottom right entry of the M-matrix is precisely the diagonal Green’s function. From (13.IL 
the time evolution of the bottom right entry is 


dtM 2 ^ 2 ix,t', z) = - 2 {dxu){x,t)M 2 , 2 {z,t) + 4:{u{x,t) + 2 z)Mi_ 2 (x,t; z). 
In particular, for z = p.j{x,t), from (13.3|) . 

(3.4) ^*^ 2 . 2 ( 2 ;, t; z) I = -2aj{x,t){u{x,t) +‘2fJ-j{x,t))- 


At any (x,t) such that pj{x,t) G (Ej ,E^), by the implicit function theorem, fj,j{x,t) is differentiable and 

dyij _ dG/dt\z=^. 
dt dG/dz\^=^,.' 

Using dSH) and (1^ . this becomes 


(3.5) 


dt 


^<Jj(u + 2pj) 



Lj) {Tj 


- Aij) n 

1^3 


(Ef 


Uj)iEl^ Hj) 

in - UjY 


1/2 


Combining this with the trace formula u = Qi o ip, we obtain an expression for dpj/dt in terms of the 
Dirichlet data. 

This shows that fij{x,t) is continuous in (£’“,£’+). By looking at the time evolution of the off-diagonal 
entries of M{x, t', Uj {x, t)), we also see that aj{x,t) is constant in t as long as fj,j{x,t) G {E~ ,E'^). 

Then, from (12.11) . (12.21) . we conclude that pj{x,t) is also differentiable in t and 


dLj 

dt 






which, when combined with (13.51) and solved for concludes the proof. 


□ 


Lemma 3.4. pj(x,t) is continuous in t. 

Proof. Due to the previous lemma, we only need to address continuity of pj at points with pj G {E~, E^}. 
Let us assume pj{xo,to) = E ^; the other case is handled analogously. 

For any e > 0, G(xo, xq; A'j'' — e,to) < 0. By (j3.1l) . the diagonal Green’s function is continuous in t; thus, 
there is a i5 > 0 such that G(xo, Xq; — e, t) < 0 for |t — to| < <5. Thus, Pj{xo, t) > E'^ — e for |t —to| < ^- CH 


Before we prove that pj is differentiable everywhere, as an intermediate step, we have to prove statements 
which will eventually be superseded. 

Lemma 3.5. Fix x G K and j G J. Lift pj(x,f) to a continuous function pj : —>■ K. 
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(i) For any a < b, 

(3.6) \ipj{x,b) - (pj{x,a)\ < / \Ej{ip{x,T))\dT. 

J a 

(a) For any a < b, if Ej{(p(x,T)) > 0 for r G {a,b), then 


(3.7) 


0 <(pjix,b) - (pj(x,a) < / Ej(ip{x,T))dT. 

J a 


Proof. We supress the parameter x throughout this proof, since it can be considered fixed. 

By uniform continuity, the interval [a, b] can be split into finitely many closed intervals I with the property 

(3.8) \Fj{u) - Fjiv)] < TT if {u,v) C I. 

By the triangle inequality, it suffices to prove (13.6L p.7l) for (a, b) C I with an interval / with the property 

( 1 ^ . 

If (pjir) ^ ttZ for r G (a, &), then p.2l) holds on (a, &), so 

Fjib) - Fj{a) = / 2j((^(T))dT, 

J a 

and dSS, (E21) follow. Otherwise, take 

ai = inf{T G {a,h)\ipj{T) G ttZ}, hi = supjr G {a,h)\ipj{T) G ttZ}. 

Then (13.21) holds on (a, oi) and (bi, b); moreover, by (13.81) and ipj{ai), (pj(bi) G ttZ, we have </5j(ai) = (pj{bi), 
so 

pai pb 

iPj (6) - ipj (a) = (fij {b) - (pj {bi) + (pj (ai) - Pj (a) = / Ej {p{T))dT + / Ej {p{T))dT 

J a J bi 

which implies (13.61) . (13.71) . □ 


Proof of Provosition For differentiability, given Lemma 13.31 there are two cases left to consider. 

First case: pj{xo,tti) G ttZ and Ej{p{xo,tQ)) = 0. Since limi-_).tQ E!j((/3(a;o,r)) = Ej{p{xo,to)) = 0, (13.61) 
implies that 

F3{xo,t) - pj{xo,to) ^ ^ 
t—^tQ t - to 

SO pj{xo,t) is differentiable at t = to and dtPj{xo,to) = 0 = Ej{p{xo,to)). 

Second case: pj{xQ,to) G ttZ and (i^(a:o, to)) ^ 0. Without loss of generality we assume Ej(p{xQ,to)) > 
0 (the other case is treated analogously). There is then a neighborhood (x, t) G (xq — e, xo + e) x (to — e, to + e) 
in which Ej{p{x,t)) > 0. We can then conclude that 


Pj{xo,to) < Pj(x,t) < pj{xo,to) + TT, if X G (xo,xo + e),t G [to, to + e); 

the first inequality follows by the strict monotonicity of the translation flow and dSH), and the second by 
joint continuity in (x, t) and shrinking e if necessary. 

For X G (xo, Xo + e) and t G (to, to + e), we can therefore use (13.21) to conclude that 

(3.9) pj{x,t) - pj{x,tQ) = Ej{p{x,T))dT 

Jto 

for X G (xo, Xo + e). Letting x —>■ xo, we conclude that (13.91) holds also for x = xo. It then follows that 


Fj{xo,t) - pj{xo,to) 
t — tf) 


^j{p{xo,to)). 


A similar argument starting with x G (xo — e,xo) and t G (to — e,to) provides the left derivative, which 
completes the proof. □ 
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4. Construction of Reflectionless Potentials and KdV Solutions by Finite Gap 

Approximants 


In this section, we construct solutions to two questions. The first result is a converse of Craig’s results, 
where we show that any Dirichlet data / G T’(S') corresponds to a reflectionless operator Hw with spectrum 
S'; in other words, we show that the map B defined by (j2.6l) is onto. In the second result, we construct 
solutions of the KdV equation with a prescribed spectrum S and prescribed initial Dirichlet data. 

On the torus 1D(S) of (angular) Dirichlet data given as before, we introduce for fc = 1, 2, 3, 

(4.1) Qk(^) = E^ + 

Lemma 4.1. If S obeys (11.161) . then Qi, Q 2 , Q 3 are continuous scalar fields on T){S). 


Proof. Since ©(S) is equipped with the product topology, it is sufficient to establish that the series in (14.11) 
converges absolutely and uniformly. For fc = 1, this follows from 

\E- +E+ - 2^j \ < 7j. 

For fc > 2, note that for x € [E~,E^], 




k 

E 

m—l 


\E-\>^-^' 


ix-E-r < 


k 

E 

m—l 


\E7 


7, 


Writing A = K + rjj^ and expanding. 


k k—r. 


m—l n—0 


k\ fk — m 
mj V n 


k-m-n n m 

Vjfilj ■ 


^,m ^ pm—1 


k — 1 — n 

utt’l. V + -jrrr^’ 


..k-1. 


k* - (E7)‘i < E E 

By the arithmetic-geometric inequality, 

where F = sup^^jy^. We therefore obtain 

k'=-(F;-)'=|<C(i7|-S,+7,) 
where C is a constant depending only on fc, S, and F. Thus, 

\iE -7 + {E+7 - 2{f,,7\< + 7 ,) 

and the series converges absolutely and uniformly by dSD. 


□ 


Our first result is a converse of Craig’s results: we show that any Dirichlet data in ©(S') corresponds to 

a reflectionless potential ©(S), or in other words, that the map B : ©(S) —>■ T>{S) defined in (12.61) is onto. 

Combined with Craig’s theory, this shows that R is a homeomorphism. 

For this result, we will assume that S is equal to its essential closure, i.e., that 

(4.2) Vx G S Ve > 0 Leb(S D (x — e, x-I-e)) > 0, 

where Leb denotes Lebesgue measure. This certainly holds if, for instance, S is the a.c. spectrum of some 
self-adjoint operator, so it will hold in our applications. 

Proposition 4.2. Let S obey (11.131) . (14.21) . Let f G ©(5) and let (f) : ^ ©('5) be the unique solution of 

the initial value problem 


dx(j){x) = 4'(()i(x)), 

m = /■ 
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If we define VF : R. —?> M 6?/ 

(4.3) VF = Qio0, 

then W G 3^(5') and B(W) = f. Moreover, if S also obeys ()1.16l) . then W G C'^(R) fl and W obeys 

the higher order trace formulas 

(4.4) Q2 0(I) = -^dlW + W^ 

Q Q IK 

(4.5) Q^ocf=-dtW--WdlW--id,Wf + W^ 

Remark. The claim that B{W) = f is not trivial; since Qi is far from injective, it is not a priori obvious 
that two different trajectories (f, (f> under the flow ip may not lead to the same potential, Qi o 4> = Qi o (j). 
The above proposition rules this out. 

The other main result of this section concerns existence of a solution of the KdV equation with prescribed 
initial Dirichlet data. 

Proposition 4.3. Let S obey (I1.15|) . (11.161) . (|4.2I) . Let f G ©(S'). There exists a function tp : —>■ ©(S) 

such that y)(0, 0) = / and 

dj;(p{x,t) = 'if{ip{x,t)), dt<f{x,t) = E{ip{x,t)). 

If we define u \ ^ M. by 

(4.6) u = QiO(p, 

then the function u{x, t) obeys the KdV equation (II.ip . Moreover, for each t gM., we have u{-, t) G ©(S) and 
B{u{-,t)) = ipfl),t). 

We will prove these propositions by a Cauchy sequence argument. The main ingredient is the theory of 
finite gap potentials, already mentioned in the introduction. 

In this and the next section, we will find it notationally convenient to assume that gaps are indexed by 
positive integers, i.e., that J = N. This doesn’t reduce generality, since in the special case where J is finite, 
the results are already known. Denote 

S^ = [^,oo)\ \J{E-,Ep 

j<N 

and denote the corresponding isospectral torus by ©(S^) and vector fields corresponding to translation and 
KdV flows by and 

Starting from an element / G T>{S), we observe the functions ip^ : —?> ©(S'^) which solve the equations 

(4.7) = = 

and obey the initial conditions 

= j<N. 

This uniquely determines the functions (p^, since the vector fields and commute, which follows 
from a direct calculation that for all i,j < N, 

Vkf - Sf 9,41^ = 0. 

At this point we recall that, by [GHl Theorem 1.48, Lemma 1.16], the function un = o ^ R. has 

the following properties: 

(i) For every t, the function UNix,t) is almost periodic in x and the Schrddinger operator has 

spectrum 

(ii) Un solves the KdV equation dm 
(hi) The Dirichlet data of uat are 
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(iv) The following trace formulas hold: 

( 4 . 8 ) Q^o^^ = un 

( 4 . 9 ) Q 2 o =-\dluN + 

Q Q 1C 

(4.10) oip^ = —dluN - -umOIun - -^{dxUNf + u%, 

lo z lo 

where are scalar fields on D(5'^) given by the expression 

N 

Q^{^) = E^ + + (E)' - 2 (^ 1 )")- 

i=i 

We begin with a statement about stability of trajectories of a vector field. 

Lemma 4.4. Assume that U, U are Lipshitz vector fields on T>{S), with Lipshitz constants less or equal to 
L. Consider solutions (j),f ®('S') of dx4> = U{(j)), dx4> = U{f). Then 

(4.11) Mx) - ^(x)lb(s) < 2 (||</)(0) - ^(O)b(s) + C) 

where m = 2L In 2 and 

C = j supyj/^ min(27r, \\Uj - Uj\\). 

E jeJ 

Proof. Let I = 1/{2L) and define the map A from C^O, Z], ©(S')) to itself by 

H9){x)=f+[ U{g{y))dy. 

Jo 

Then ^ 

||A( 5 )(a;) - A(/i)(a;)|| < f \\U{giy)) - U{h{y))\\dy < Lx\\g - h\\ 

Jo 

so, by our choice of /, 

||A(g)-A(/i)||<i||g-/i||. 

Since A is a contraction with coefficient 1/2, it has a unique fixed point, which is precisely the path 4>{x). 
Similarly, defining A on (^([O, Z], ©(S)) by 

H9){x) = f+[ U{g{y))dy, 

Jo 

A is a contraction with coefficient 1/2 whose unique fixed point is (j){x). 

There are two ways of estimating Aj{g) — Aj{g): one is by integrating, 

\\U9)-U9)\\<\\fj-m+m-Uj\\, 

and the other is the trivial estimate ||Aj(g)—Aj( 5 )|| < tt. Using the smaller of these two estimates, multiplying 
by and taking sup^gj, we conclude 

||A(g)-A(5)||<||/-/|| + ^. 

In particular, for g = 4>, 

By standard estimates from fixed point theorems, since f = lim„_,.oo A"(^) and A is a 1/2-contraction, this 
implies that 

||/>-0||<2||/-/||+C, 

that is, 

Mix) - f{x)\\ < 2\\(j){0) - ^(0)11 +C, Vx G [0,Z]. 
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Applying this inductively, we see that for n G N, 

U{nl) - knl)\\ < 2" (ll^(O) - .^(0)11 + c) - C 

and therefore for a; > 0, 

||<^(cr) - ,^(a:)|| < (||</>(0) - ^(0)|| +c)-C. 

The same argument works for negative x by reversing the vector fields, so we conclude (14.111) . 


□ 


To apply the previous lemma to compare trajectories on different isospectral tori, we will lift all the fields 
and solutions to T){S) in the following way. 

There is a natural projection tt : Ti{S) —>■ given by 

= (t>j, j < N. 

To lift the scalar fields, simply define : T){S) —>■ K by 
Introduce the vector field 4'^ on D{S) by 


Qk — Qk ° 




J<N 

3>N 


and, analogously, define the vector field 2^ on ©(S') by 


( 


j<N 


(7r((^)) j > N 


To define —>■ ©(S), set the initial value 

<^^( 0 , 0 )=/, 

determine its values for x = 0 by requiring 

t) = t)), 'it e R, 

and then its values for arbitrary x by requiring 

{x,t) — {x,t)), Vx,tsR. 

With these definitions, obviously, 

N ~N 

If = TT O ip . 

We similarly introduce 93 : R^ —^ ®(S) by 

<^( 0 , 0 ) = /, 

dt(p{0, t) = 2((/3(0, t)), it e R, 
dx<p{x,t) = it{(p{x,t)), Vx,tsR. 

Lemma 4.5. Assume that the set S obeys (11.151) . There exists m > 0 and constants Kjq such that 
limAr->.oo A'at = 0 and, for all x, t G R, and all j < N, 

(4.12) \\T^ix,t) - pix,t)\\ < 


Proof of Lemma \4-.5\ If (jl.l5p holds, we know that ik, 2 are Lipshitz vector fields; the proof of these facts 
gives explicit upper bounds for the Lipshitz constants in terms of gap sizes and distances, so it also applies 
to and 2^, giving uniform Lipshitz estimates in N. Denote by L such an upper bound on the Lipshitz 
constant which works for ik, 2, and all the tk^, 2^. 

In particular, by the Lipshitz property, values of p>{x,t) and (p^ (x,t) are uniquely determined by the 
above definitions. 
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For j < N, 


1]/ - _ 


< l^f I 


n 


i-^k Mi) 


(Mfc - Mi)^ 


1/2 


- 1 


< 2(?7i.o + 7i)^^^ n 


l<fc<A 

fc#i 


Ik 

Vj,k 


1/2 


n 1+ 


^fe=IV+l 


Ik 

Vj,k 


1/2 


- 1 


< 2{Cj - Cj^n) 


where Cj is defined in (jl.l2l) and 


Cj,N — {r]j,o +lj) n ( 1 + ~ 

1<N ^ 


1/2 


Since = —2{Qi + 2/ij)^'^, we estimate 


■3J^j 
^ 

“i “i 


where 


Let m = 2L In 2 and let 


< 2|Qi - Qf I I + 2|Qi + 2/1,1 Id/, - d/, 

< {Dj - Dj^N)Cj^N + Dj{Cj - Cj^n) 

^ DjCj — Dj^j^Cj^jsf 


N 

L>,-,Ar = 2V'7i +4?7,_o + 4|S|, Dj = lim 

• ^ M —Vf-yn 


N\ 


Z=1 


Kn = j- sup max (27r, Cj — C/.at, DjCj — Dj^NCj^N) ■ 


j<N 


Since lim,_>oo 7/ = 0 and, for each j, Cj^N —>■ Cj, Dj^jq Dj as iV —i oo, it is clear that limAr->.oo Kn = 0. 
Applying Lemma [4.4l to compare trajectories of vector fields and S with initial condition /, we conclude 

(4.13) - ^(0,t)|b(s) < 

Applying Lemma to compare trajectories of d'^ and d/ using (14.131) as initial conditions, we obtain 

Mx,t) - g3^{x,t)U^s) < 2(KAe™l*l + 

which implies (I4.12L □ 


If we remove all discussion of t-dependence and of the vector fields S, from the previous proof, then 
we only need the Lipshitz property to hold for d/, d/^, and the argument only requires the weaker condition 
(11.131) instead of (|1.15|) . Likewise, the constant Kn from the previous proof can be redefined as 

O 

Kn = T sup 7 ^^ min(27r, Cj - Cj n)- 

C j<N 

The argument then yields the following result. 

Lemma 4.6. Assume that the set S obeys (11.131) . There exists m > 0 and constants Kn such that 
lim^v-j-oo Kn = 0 and, for all a: G K., and all j < N, 

(4.14) 0) - (fi{x, 0)11 < 
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Proof of Prop, By Lemma 321 

P(x,0) = lim 

N—^oo 

converges uniformly on compact sets of x. Denote 

W{x) = Qi(vj(x,0)). 

Since are uniformly continuous in N, we conclude that the sequence {x) := ui^{x,0) converges 
uniformly on compacts to W {x). 

We will now use a technique of Craig m Section 5]; in [Cr] it is applied to a limit, uniform on compacts, 
of a sequence of translates of a given potential, but it applies equally well to the limit of the sequence . 
Since the Schrddinger operators HyyN converge in strong resolvent sense to Hw, and since a{H^^N) = , 

we conclude that a(Plw) C S. Moreover, uniform convergence on compacts also implies convergence of 
the diagonal Green’s function uniformly on compacts (in z) away from S. Since the are reflectionless, 
fiG{x, x-,z -\- jO, W^) = 0 for a.e. z G S', so by [Crl Lemma 5.2], lkG{x, x; z + *0, W) = 0 for a.e. z € S. 

By analyticity, SG(a;, x; z, W) = 0 for z G M \ a{P[w)- Since G{x, x; z + iO, W) = 0 is possible only on a 
set of Lebesgue measure 0, combining this with the previous conclusion shows that S \ a(Plw) has Lebesgue 
measure 0. Since S and a{Hw) are closed sets, (14.21) then implies that <7{P[w) = S and W G 3?(S). 

We will now prove that B{W) = /. Denote by ^Ar(a:), (x) the Dirichlet data for (related by (12.11) . 

(12.21) to the angular Dirichlet data ip^{x)), and by ^(t), cr{x) their limits. Fix j G J and recall that, by the 
way were constructed, for N > j, G{x,x; z,W^) > 0 for z G [iJ,^(x),E^) and G{x,x-, z,W^) < 0 for 
z G {E~, (x)]. Using again convergence of the diagonal Green’s function, we conclude that G{x, x; z, W) > 
0 for z G [pLj{x),E'^) and G(x,x;z, lU) < 0 for z G [EJ, pj{x)\, so ^j(x) are precisely the Dirichlet data 
of W. Since crj{x) can be read off from whether /ij (x) is increasing or decreasing, we conclude that the 
Dirichlet data of W{x) are precisely <f{x). In particular, B{W) = i^(0) = /. 

If (|1.16l) also hold, then are uniformly continuous, so Q 2 o and o are uniformly 

Gauchy on compacts. 

Fix a compact K C M.. Since we already know that the sequence is Gauchy in C{K), from (14.91) we 
conclude that df.W^ is also Gauchy in G{K), so is Gauchy in G^{K). Now, similarly, from (I4.10|) we 
conclude that d'^W^ is Gauchy in C{K), so is Gauchy in C^{K). Its limit is W, so we conclude that 
W G C{K) and, by taking limits of (14.91) . (I4.10p . that (14.41) . (14.51) hold. 

Finally, by using boundedness of Qi,Q 2 ,Q 3 , we conclude from (14.31) . (14.4L (14.51) that W G IU‘*’°°(R). □ 

Proof of Prop, \4-3l By Lemma |4.51 convergence 

(4.15) (p{x,t) = lim (f^{x,t) 

N^oo 

is uniform on compacts in (x,t) G Introduce 

(4.16) u{x,t) = Qi{(p{x,t)). 

By Prop. 1121 (I4.I6|) implies that u{-,t) G fk(5') for each t and B{u{-,t)) = (p{0,t). 

In Prop. 14.21 it was proved that for each t G R and each compact iF C R, 

(4.17) lim |luiv(-,i) - u(-, t)||c4(A) = 0 

N—yco 

In the setting of this proposition, we know that convergence (14.151) is uniform on compacts in R^, so conver¬ 
gence (14.171) is uniform on compacts in t. In particular, this implies that d^u is jointly continuous in (x, t) 
for fc = 0,1,2,3,4. 

Since un obey the KdV equation dni), 

ziAr(x, t) = WAr(x, 0) + / {6uNdxUN — dx'^N)ix,r)dr. 

Jo 
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By uniform convergence, we can take the limit N ^ oo to conclude 

u{x,t) = u{x,0) + / {6udxU — d^u){x,T)dT. 

Jo 

Since the integrand is continuous, this implies that u{x,t) is differentiable in t and obeys (11.11) . □ 


5. Linearization of the KdV Flow 


In this section, we will resume our investigation of the finite zone solutions U]si{x,t) of the KdV equation 
and their limit u{x, t). We will consider how these behave with respect to the generalized Abel map introduced 
by Sodin-Yuditskii. We assume throughout this section that S' is a regular Parreau-Widom set with finite 
gap length. 

In order to construct a reparametrization of the isospectral torus, Sodin-Yuditskii [SYl] introduced func¬ 
tions ^j{z), one for each open gap {E~ ,E^ ). The map Cj{z) is the solution of the Dirichlet problem on C \ S 
with boundary conditions on S given by 


(5.1) 


ijix) 


1 a; = ooora;GS,a;> E^ 
0 X € S,x < EJ 


Regularity of S implies that (z) is a continuous function on C and a harmonic function on C \ S, and its 
values on S are given by dsn). 

Let 7 r(C \ S) be the fundamental group of C \ S. It is a free group with the set of generators given by 
{cj}jgj, where Cj is a counterclockwise simple loop which intersects R. at the points A — 1 and {E~ + E^)/2. 

Following [SYlj . consider the group 7 r*(C\S) of unimodular characters of 7 r(C\S'). Use additive notation 
for 7 r*(C \ S). An element a € tt* {<C \ S) is uniquely determined by its action on loops Cj, so we can write 
a = {aj}j^j where aj = a{cj) £ T. Endow 7 r*(C \ S) with the topology dual to the discrete topology of 
7 r(C \ 5"); there are many ways to choose a metric which induces this topology, for instance 


d{a, a) = min(|aj — 5^ |, 7 ^), a,a G tt*{C \ S), 

3&J 


where 7 ^ = E^ — Ej , as before. 

[SYl] define the Abel map A : T>{S) —>■ 7 r*(C \ S) by defining its components Aj = A(cj), 


(5.2) Aj{(p) = TT^cTfc (O(Mfe) - )) (mod 27rZ) 

keJ 

where, as always, we assume that /i^, (Tk are given in terms of ipk by (lO) . (IQ) . 

[SYl] prove the following (assuming that S' is a regular Parreau-Widom set with finite gap length): 

(i) For each j, the sum in (15.21) converges absolutely and uniformly in £ ®(5'); in particular, the map A 
is well-defined; 

(ii) A is a homeomorphism between D(S) and 7r*(C \ S); 

(hi) Denote by A^ : ©(S^) —>■ 7r*(C \ S^) the Abel map for . Project R(S) to D(S^) by truncation, 

and embed 7r*(C \ S^) into 7r*(C \ S) by assuming (cj) = 0 for j > N. With those conventions, 
consider A^ as a map from 'D(S) to 7r*(C \ S). Then, A^ —>■ A as iV — 00 , uniformly on T)(S). 
mu introduced this map to generalize the notion of Jacobi inversion, which exists in the finite-gap 
setting. In that setting, it is known (see, e.g., [GHl Thm. 1.44]) that Jacobi inversion linearizes translation 
and KdV flows. Therefore 


(5.3) A^i^^ix, t)) = A^i^^iO, 0)) + S^x + Ct 

for some ^ ■ 

We can now define the map 

M = B-^o : 7r*(C \ 5') ^ 31(5'). 
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Proposition 5.1. Assume that S is a regular Parreau-Widom set which obeys (11.151) . (11.161) . (14.2p . Then, 
the map M, is a homeomorphism if'Ji{S) is equipped with the metric inherited from VF'^’“(]R). 

Recall that by Prop. 1121 we already know that Ik(S') C 


Proof. For any ei > 0, there exist e 2 , £3 >0 such that 

sup||(a + Sx) - (a + fe)||^.(c\s) < £3 

xGR 

sup||A“^(a + Sx) - A~^{a + fe)||D(s) < £2 
sup sup|Qfc(A“^(a + Sx)) - Qk{A~^{a + fe))| < £1 

feG{1.2,3} xGR 

Denote W = M(a), W = M(a). By (14.311 . the above estimate for fc = 1 implies 

IIVF-T^lloo <£i, 

so M is continuous as a map to L°° (R). Inductively, using (14.41) and (14.51) , we conclude that M is a continuous 
as a map into IF^’°°(R) and IF'^’°°(R), respectively. □ 


\a - a||,r*(c\S) < £3 


Proof of Theorem[M Since Pfy is almost periodic and a{Hv) = CTac(^Iv) = S', a result of Remling [Rel] 
implies V G 3?(S). Therefore, V corresponds to a set of Dirichlet data / = B{V), as defined in (12.6p . 

By results of Rybkin |Ry| reviewed in Section 3, if u{x,t) is a solution of (11.11) . (11.21) which obeys (II.7L 
then u{-,t) G lk(S) for each t G R, and the Dirichlet data (p{x,t) of u obey 

(5.4) d„ip = '^{ip), dtip = E{ip) 


and (/j(0, 0) = /. This determines (p{x,t) uniquely, since 111,2 are Lipshitz vector fields. The trace formula 
u = Qi o If then determines u uniquely. 

Conversely, by Prop. 1231 the solution (p{x,t) of (15.41) with V5(0,0) = / exists and generates a solution of 
(11.11) hy u = Qi o ip. Since R, w(-,0) G 1R(S) and B{V) = R(u(-,0)) = /, we conclude that V = u(-,0), so u 
obeys the correct initial condition m- We have therefore established existence and uniqueness. 

Since Hy has purely a.c. spectrum, the set S obeys dOl), so Prop. [5Tl applies to S. 

We now recall the functions (p^{x,t) introduced in Section 4. Since (p^ —> (p uniformly on compacts 
and A^ —?> A uniformly, we can conclude that A^{ip^{x,t)) converge uniformly on compacts to A{ip{x,t)). 
Taking the j-th component of (15.31) . it follows from uniform convergence that the limits 


Sj = lim Si, 

•' N^oo •’ 


N 


0 = lim C 


N 


N —¥oo 


exist and 

Aj{ip{x,t)) = Aj{ip{Q, 0)) + SjX + Qt. 
In particular, = M“^(R) + C,t. 


□ 


6. Application to Small Quasi-Periodic Initial Data 

In this section, we will show that Theorem [1] follows from Theorem [31 This will rely on spectral properties 
of Schrodinger operators with small quasiperiodic initial data, which have been extensively studied in [Ell 
lDCTTllDnblllD(IF2irDCTn3] . 

Until now, we have indexed gaps by the abstract gap label J, and in some sections it was notationally 
convenient to assume J to be equal to N. However, for almost periodic Schrodinger operators, there is a 
natural gap label given by the rotation number, see [JoMoj . In the setting of Theorem [I] this means that 
gaps are naturally labelled by m G with m and —m corresponding to the same gap and with m = 0 
corresponding to the bottom of the spectrum. We will use this label from now on. We remind the reader 
that in this paper, for m = (mi,..., m^,) G IT, \m\ stands for the norm, |m| = X]j=i|mj|. 

In the setting of Theorem[Tl it is known that there is an £o(ao, 5o, kq) > 0 such that, if £ < £ 0 , to obeys 
the Diophantine condition and V G lP(w, £, kq), then, with S = a{Hv): 
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(i) Hy has purely a.c. spectrum; 

(ii) For every m € \ {0}, 

(6.1) 7 m < 2eexp (“ 

(iii) For every m € 1/ \ {0} and n £1/ with m ^ n and \m\ > |n|, we have 

(6.2) ?7m.n > 

for some constants a,b > 0 depending only on oj, e, kq, v. 

(iv) For every m G \ {0}, 

(6.3) T]m,o < c|TOp, 

for some constant c > 0 depending only on w, e. 

(v) The set S is homogeneous in the sense of Carleson: (11.111) holds for r = 1/2. 

(vi) If W G (K(S'), then W is also quasiperiodic with the same frequency uj and similar decay properties of 
Fourier coefficients; more precisely, W G lP(w, \/4e, ko/4). 

We will need a subexponential estimate of the constants 

Cm = (7m.0+7m)'/" n + ^ ' 

n^m 

Fix a constant L = L(e, kq, a, b) such that 

2eexp Vm G Z'^ \ {0}. 

Define 

i?m = {n G Z"" \ {0} I n yf to, 7 „ > 

It follows from (|6.1F (16.21) . and our choice of L, that n G Pm implies \n\ < \m\. In particular, the set Rm is 
finite. We order its elements by decreasing norm, i.e., we denote Rm = with 

(6.4) |n®| > ••• > 

Lemma 6.1. (i) There are constants t, /3 which depend only on a,b,KQ,v, such that 

(6.5) I > T exp |) , Vr < tq - 2. 

(ii) There is a constant D = D(a, 6, kq, v) such that 

\Rm\ < log 2 log 2 |TO| + D. 

(iii) There is a constant F = F(a, 6, kq, u, w) such that 

(6.6) Cm < F’exp(Flog|TO| loglog|TO|). 

Proof, (i) Among rim^, nm"*"^^, nm"*"^^, we can find two for which the corresponding gaps lie on the same side of 
the TO-th gap; denote them s = nm^\t = n^\ with ri < r 2 . By the ordering (16.41) . since r < ri < r 2 < r + 2, 
it suffices to prove 

(6.7) |s| > rexp(/3|t|). 

Since the gaps labelled by s, t are on the same side of the m-th gap, 

Vs,t < max{r]m,t,V-m,s) < < L“^/‘‘(2e)^/^ exp (-ylil) 

Combining this with (16.2p gives 

a|s|"'’ < exp (-y |i|) , 

which can be rewritten in the form (I6J1) . 
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(ii) Pick R = R{a,b,KQ,i') > 2 such that Texp(/?x) > for x > R. Assume that \n'm\ > R for all 
r < 2ri. Then (|6.5I) gives 

\m\ > I > |^ > • • • > V"' > 2"^"', 

which implies 2ri < log 2 log 2 |TO|, bounding the number of with \ri^\ > R. Trivially, the number of 
with \n'm\ < R is less than (2i? + 1)^. Combining the two estimates completes the proof. 

(iii) To estimate Cm, we estimate first the part of the product with n ^ Rm, 


^ n^J\R^ VneJ 


1 < C 


n I — 


nGJ\Rm ^ ' ' nGJ\Rrr 

by (jb.ip . For the part of the product with n G Rm, use \n\ < \m\ and the bound on \Rm\ to conclude that 


tiGR 


Finally, the factor r]m,o + 7 m is polynomially bounded by (IQ) and (ICT . 


□ 


Proof of Theorem QJ In order to apply Theorem [2l it only remains to verify that Hy obeys the conditions 

(nm), (inni. 

(11.161) follows immediately from 

( 1 + ?ym,o) 7 m < 2 e(l+ a^) ^ exp (^-y |to|) < oo. 

mGZ'^\{0} meZ'^\{0} 

1 / *? 

Similarly, ^ 7 m < oo follows immediately from (|6.1I1 . The second inequality in (jl.l5|) , 

< oo. 


1/2 1 + Vm,0 ^ 

sup 7m -W 

m^J Vm,0 


holds since the exponential decay of 7 m controls the (at most subexponential) growth of 77 ^and of Cm- 
For the third inequality in (11.151) . use 


1/2 1/2 

2 =^ < {2eL)'/VJ'yP 

Vm,n 


for all m, n G J to conclude that 


sup ^ 


m^J 


n^J 

n^m 


' 1/2 1 / 2 ' 
7m In 

^m,n 


(1 + Vm,o)Cm < sup V (qmSy"^) (1 + VTn,o)Cm < OO, 

mGJ - ^ ' 


n^J 

n^m 


again by the observation that the exponential decay of 7 m controls the at most subexponential growth of 
the other factors. 

Therefore, Theoremapplies to the initial data V G lP(a;, e, kq)- This implies Theorem [H since Ik(S') C 
V{uj, \/4i, ko/4) by |DCL2| . □ 
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